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Abstract 

The Hartree-Fock (HF) approximation with the charge-density-wave (CDW) 
instabiUty is used to study the ground-state phase diagram of the spinless 
Fahcov-Kimbah model (FKM) extended by /-/ hopping in two and three di- 
mensions. It is shown that the HF solutions with the CDW instability repro- 
duce perfectly the two-dimensional intermediate coupling phase diagram of the 
FKM model with /-/ hopping calculated recently by constrained path Monte 
Carlo (CPMC) method. Using this fact we have extended our HF study on 
cases that have been not described by CPMC, and namely, (i) the case of small 
values of /-electron hopping integrals, (ii) the case of weak Coulomb interac- 
tions and (iii) the three-dimensional case. We have found that ferroelectricity 
remains robust with respect to the reducing strength of coupling (/-electron 
hopping) as well as with respect to the increasing dimension of the system. 

PACS nrs.: 71.27.+a, 71.28.+d, 71.30.+h 
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The Falicov- Kimball model (FKM) is a paradigmatic example of simple model to 
study correlation effects in interacting fermion systems on a lattice [1]. The model 
was originally proposed to describe metal-insulator transitions and has since been 
investigated in connection with a variety of problems such as binary alloys [2], the 
formation of ionic crystals [3], and ordering in mixed- valence systems [4]. In the last 
few years the FKM was extensively studied in connection with the exciting idea of 
electronic ferroelectricity [5, 6, 7, 8, 9, 10]. The motivation for these studies comes 
from the pioneering work of Portengen at al. [11, 12] who studied the FKM with a 
k-dependent hybridization in the Hartree-Fock (HF) approximation and found that 
the Coulomb interaction U between the itinerant d-electrons and the localized /- 
electrons gives rise a non-vanishing excitonic (/"'"(i)-expectation value even in the 
limit of vanishing hybridization V ^ 0. As an applied (optical) electrical field pro- 
vides for excitations between d- and f-states and thus for a polarization expectation 
value Pfd = {f^di), the finding of a spontaneous Pfd (without hybridization or elec- 
tric field) has been interpreted as evidence for electronic ferroelectricity. This result 
stimulated further theoretical studies of the model. Analytical calculations within 
well controlled approximation (for U small) performed by Czycholl [5] in infinite di- 
mensions did not confirm the existence of electronic ferroelectricity. In contrast to 
results obtained by Portengen et al. [11, 12] he found that the FKM in the symmetric 
case {uf = Ud = 0.5) does not allow for a ferroelectric ground state with a spon- 
taneous polarization, i.e., there is no nonvanishing (/+(i)-expectation value in the 
limit of vanishing hybridization. The same conclusion has been also obtained inde- 
pendently by extrapolation of small-cluster exact-diagonalization and density matrix 
renormalization group (DMRG) calculations in the one dimension for both interme- 
diate and strong interactions [6] . In these regions the finite-size effects are neghgible 
and thus the results can be satisfactory extrapolated to the thermodynamic limit. 

Hybridization between the itinerant d and localized / states, however, is not the 
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only way to develop d-f coherence. Recent theoretical works of Batista et al. [8, 9] 
showed that the ground state with a spontaneous electric polarization can also be 
induced by /-/ hopping for dimensions D > 1. In the strong coupling limit this result 
has been proven by mapping the extended FKM into the xxz spin 1/2 model with 
a magnetic field along the 2;-direction, while in the intermediate coupling regime 
the ferroelectric state has been identified numerically by constrained path Monte 
Carlo (CPMC) technique. On the base of these results the authors postulated the 
following conditions that favor the formation of the electronically driven ferroelectric 
state: (a) The system must be in a mixed-valence regime and the two bands involved 
must have different parity, (b) It is best, though not necessary, if both bands have 
similar bandwidths. (c) A local Coulomb repulsion (U) between the different orbitals 
is required. 

In the present paper we study the extended FKM (the spinless FKM with /-/ 
hopping) in the HF approximation with the charge-density- wave (CDW) instability. 
For reasons mentioned above we restrict our studies on dimensions D > 1. First 
we show that the HF solutions with the CDW instabihty reproduce perfectly the 
ground-state phase diagram obtained by the CPMC method (D = 2) for intermediate 
Coulomb interactions [9]. This "calibration" allows us to extend calculations to the 
case of small values of the /-electron hopping integral \tf\ < 0.1, that has been 
omitted in the CPMC phase diagram for numerical problems. Just in this region 
we have found a new phase that corresponds to the inhomogeneous solution for 
(/+d)-expectation value. This result completes the ground-state phase diagram of 
the two-dimensional FKM extended by /-/ hopping for intermediate couplings. The 
same calculations we have performed also in the weak coupling limit (for D = 2) as 
well as in three dimensions. We have found that the ferroelectricity remains robust 
with respect to the reducing strength of the coupling as well as with respect to the 
increasing dimension of the system. 
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1 The model 



The extended FKM for the spinless fermions on a D-dimensional hypercubic lattice 
is 

h^-uy: did, -tfY: ftfj + uy: ^t^^dtd^ + EfY: Ufi, (1) 

<ij> <ij> i i 

where fj^ (fi) and df (di) is the creation (annihilation) operator of heavy (/) and 
light {d) electron at lattice site i. 

The first two terms of (1) are the kinetic energies corresponding to quantum- 
mechanical hopping of d and / electrons between the nearest neighbor sites i and 
j with hopping probabilities and tf, respectively. The third term represents the 
on-site Coulomb interaction between the d electrons with density Ud — dfdi 
and the / electrons with density n/ = jJ2iftfi, where L is the number of lattice 
sites. Usually, the hopping integral of the d electrons is taken to be the unit of energy 
{td — 1) and the /-electron hopping integral is considered in the limit \tf\ < 1. This 
is a reason why the d electrons are called light and the / electrons heavy. 

In our HF study of the extended FKM we go beyond the usual HF approach [13] 
in which only homogeneous solutions are postulated. In accordance with [14] we 
consider here also inhomogeneous solutions modeled by a periodic modulation of the 
order parameters: 

(n/) = n^ + (5/cos(Q-r,) , (2) 
(nf)=n'^ + 5rfCos(QT,) , (3) 
= A + Apcos(Qt,) . (4) 

where Sa and 5/ is the order parameter of the CDW state for the d- and /-electrons 
and A is the excitonic average. The nesting vector Q = (tt, tt) for D = 2 and 
Q = (tt, tt, tt) for D = 3. 

Using expressions for {n{),{nf) and {fi~di) the HF Hamiltonian of the extended 
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FKM can be written as 



- taj:dtdj-tfJ2fUj + Efj:n{ + Uj2{nf + SfCOs{Q-r,))nf 

{hi) {hi) i i 

+ UY.{n^ + (5dCos(Q • v,))ni - UY,{A + Apcos(Q • r,))d+/, + h.c. (5) 

i i 

Following the work of Brydon et al. [14] the effective HF Hamiltonian is diagonalized 
by canonical transformation 

7^ = Kdk + v^dk+Q + afh + b^fk+Q , 1,2,3,4, (6) 

where a'^,bf,u^,v'^ are solutions of the associated Bogohubov-de Gennes (BdG) 
eigenequations: 



k ) 



(7) 



with 



U5f 
-UA 
-UAp 



US 



f 



-k+Q 

-UAp 
-UA 



-UA 

+ Unf -UAp 

ei + Un'^ + Ef 
USd 



and 



-UAp \ 
-UA 

U5d 

eUg + Un^ + Ef J 



(8) 



(9) 



V J 

The corresponding energy dispersions ef and e{ can be obtained directly by the 
Fourier transform of the d- and /-electron hopping amplitudes and for the case of 
hypercubic lattice they are given by {a — d, f): 



f"' — 



f"' — 



-2to,{cos{kx) + cos(fcy)), for Y)—2, 
-2ta{cos{kx) + cos{ky) + cos{kz)), for D=3. 



(10) 
(11) 



The HF parameters Ud, Sd, rif, 5f, A, Ap can be written directly in terms of the 
Bogohubov-de Gennes eigenvectors: 

- ^ E ' E{« + <vnf{ET) ■ (12) 

k m 

= ^ E'E{« + «}/(^r) • (13) 

k m 

^^-^E'E{«+ w}/(^r) • (14) 

'^Z = ^ E ' E{^r< + <b]:}fiEr) . (15) 

k m 

^ = 4 E ' E{« + bTvnnm ■ m 

k m 

= ^ E ' EiKuT + <vnf{K) ■ (17) 

k rn 

where the prime denotes summation over half the Brillouin zone and f{E) = 1/{1 + 
exp [/?(£' — //)]} is the Fermi distribution function. 

The same approach has been used recently by Brydon et al. [14] to study the 
interplay between excitonic effects and the CDW instability in the FKM with on-site 
as well as non-local hybridization. Here we use the zero temperature variant of this 
procedure to describe ground-state phase diagram of the spinless FKM with /-/ 
hopping. 

2 Results and Discussion 

To determine the ground-state phase diagram of the extended FKM in the Ef-tf 
plane (corresponding to selected U) the HF equations are solved self-consistently for 
each pair of {Ef,tf) values. We use an exact diagonalization method to solve the 
Bogoliubov-de Gennes equation. We start with an initial set of order parameters. 
By solving Eq. (7), the new order parameters are computed via Eqs. (12) to (17) and 
are substituted back into Eq. (7). The iteration is repeated until a desired accuracy 
is achieved. 
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First we have examined the two-dimensional extended FKM model in the inter- 
mediate coupling regime and tf negative. For this case there exists the comprehensive 
phase diagram of the model obtained by a CPMC technique [9] for /-electron hop- 
ping integrals \tf\ > 0.1. According these Monte-Carlo studies the phase diagram of 
the extended FKM consists of only three main phases, and namely, (i) the integer- 
valent state {uf — 0, 1, rid — 1, 0), (ii) the mixed-valent CDW state (n/ — na — 0.5), 
and (iii) the mixed-valent ferroelectric state that is stable for remaining values of n/ 

In Fig. 1 we have displayed typical examples of our HF solutions obtained for 
n(i,Sd,nf,Sf, A, Ap in the intermediate coupling regime U — 2. It is seen that the 
extended FKM in the HF approximation with the CDW instability exhibits non- 
vanishing excitonic (/"'"(i)-expectation value for all /-electron densities except the 
case when Uf — 0, 1/2 and 1. Thus in accordance with the Quantum Monte-Carlo 
studies [9] we have found that the ferroelectric ground state with the spontaneous 
polarization is stabilized when the system is in the mixed valence regime and the sign 
of the /-electron hopping integral is opposite to the sign of the d-electron one. The 
fact that HF solutions can describe the existence of ferroelectric ground-state with 
spontaneous polarization is not surprising, since this state has been found already 
in the homogeneous HF solution of the conventional FKM {tf — 0) in the limit of 
vanishing hybridization V [11, 12], even for all /-electron concentrations (for 
all values of from the d-electron band). What is however surprising is that the 
HF solutions with the CDW instabihty reproduce perfectly the ground-state phase 
diagram obtained by CPMC method for all examined values of /-electron hopping 
{\tf \ > 0.1). This is clearly demonstrated in Fig. 2, where both phase diagrams are 
compared. 

The fact that the HF approximation with the CDW instability can describe quali- 
tatively as well as quantitatively ground-state properties of the FKM with /-electron 

7 



hopping motivated us to extend our HF study on cases that have been not described 
by Quantum Monte-Carlo simulations. At first this is the case of small /-electron 
hopping integrals {\tf\ < 0.1) that has been not considered in the original work 
of Batista et al. [9] because numerical difficulties which appear in the Quantum 
Monte-Carlo simulations for small tf (the limitations in the numerical accuracy). 
The second interesting case that we would like to study here within the HF theory 
is the three-dimensional case for which the numerical results are very rare due to 
numerical limitations on the size of clusters. 

Let us first discuss our two-dimensional results obtained in the limit of small 
values of /-electron hopping integrals. In Fig. 3 we present results of detailed HF 
analysis performed in this limit for A, Ap and dd- It is seen that the non- vanishing 
excitonic {f~^d) expectation value persists also for small values of \tf\ but now the 
inhomogeneous solution Ap (with AB-sublattice oscillations in the excitonic 
and charge order parameters) is stabilized against the homogeneous one (Ap = 0). 
The effect is especially strong when we approach the tf = limit. This is clearly 
demonstrated in Fig. 4, where the complete intermediate-coupling phase diagram of 
the FKM with /-/-hopping is displayed. Five different phases depicted in Fig. 4 as 
a (the full / band), (the excitonic phases), 7 (the CDW phase) and e (the full 
d band) correspond to following HF solutions: 



a phase: 


A = 0, 


Ap 


= 0, 




= 0, 




= 0, 


Uf = 1 






P phase: 


A > 0, 


Ap 


= 0, 


Sf 


= 0, 


Sd 


= 0, 


< Uf < n'j 


for Ef 


> 


















l-nj <nf <1 


for Ef 
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(3' phase: 
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<0, 
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Sd 
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n'j < Uf < 1/2 


for Ef 


> 


















1/2 < Uf < 1 - nj 


for Ef 
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7 phase: 


A = 0, 


Ap 


= 0, 


5/ 


<0, 


Sd 


>0, 


Uf = 1/2 






e phase: 


A = 0, 


Ap 


-0, 




= 0, 


Sd 


= 0, 









The stability of different HF solutions was also checked numerically by calculating 
the total energy and it was found that all phases presented in the ground-state phase 
diagram represent the most stable HF solutions. To determine the type of transitions 
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between different phases we have performed an exhaustive numerical study of the Ef 
dependence of the HF order parameters (the typical examples are shown in Fig. 1 
and Fig. 3). At a first glance it seems that there are both first-order {tf large) and 
second-order {tf small) phase transitions in the extended FKM with /-/ hopping. 
However, a more detailed analysis of numerical data (with much higher resolution 
than used in Fig. 1 and Fig. 3) showed that the /3' phase persists also for large tf, 
although its stability region is now considerably reduced (see insets in Fig. 1). Thus 
there is no difference between the case of small and large values oitf. In both cases 
the HF order parameters change continuously indicating that the phase transitions 
between different phases presented in the (Ef-tf) ground-state phase diagram are of 
the second order. 

The same calculations we have performed also in the weak coupling limit {U < 1). 
We have found that the phase diagrams obtained in the weak and intermediate 
coupling regime have qualitatively the same form and only one difference between 
them is that the ferroelectric domain (/3) is stabilized against remaining phases with 
decreasing Coulomb interaction (see Fig. 4). Of course, this fact does not imply 
automatically that the excitonic (/^rf) expectation value persists also for vanishing 
U and that the Coulomb interaction U is not necessary for a stabilization of the 
ferroelectric state, what should be in a contradiction with conclusions based on the 
CPMC simulations. Indeed, calculations that we have performed for different values 
oitf at the selected /-electron density n/ = 1/4 showed (see Fig. 5) that the excitonic 
(f'^d) expectation value is zero for U — 0, rapidly increases with increasing U and 
tends to the saturated state for U sufficiently large. This confirms independently the 
third postulate of Batista et al. [9] and namely, that the local Coulomb interaction 
between the different orbitals is required in order to stabilize the ferroelectric state 
with the spontaneous polarization. 

Before discussing the case of positive tf let us show explicitly the HF solution 
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for the limit of the conventional FKM {tf — 0). For this case we have found that 
A = Ap = in the 0,7 and e phase, while A = — Ap in the phase. The last 
solution implies that the extitonic (/j^(ij)-expectation value is equal to 2A on the 
A sublattice of the hypercubic lattice, while {f^di) = on the B sublattice. For 
the symmetric case Ef = our solutions are fully consistent with the CzychoU's 
ones obtained in the limit of infinite dimensions [5]. On the other hand both these 
inhomogeneous solutions fully differ from the homogeneous one [11, 12] that predicts 
a non-zero excitonic (/j^cij)-cxpcctation value for all Ef from the mixed valence 
regime with maximum of {fi~di) at Ef — 0. 

Similar calculations as for < we have performed also for i/ > 0. We have 
found that the ground-state phase diagram for t/ > has exactly the same form as 
for tf < 0, however five different phases a, (5, 7 and s are now characterized by: 



a phase: 


A = 0, 


Ap 


= 0, 




= 0, 


Sd 


= 0, 


Uf = 1 






P phase: 


A = 0, 
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= 0, 
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= 0, 
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< 


phase: 
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1/2 <nf <l-rff 


for Ef 
for Ef 
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7 phase: 


A = 0, 


Ap 


= 0, 




<0, 


Sd 


>0, 


Uf = 1/2 






e phase: 


A = 0, 


Ap 


= 0, 




= 0, 


Sd 


= 


n/ = 







Thus the main difference between the phase diagrams obtained for negative and 
positive tf is that the ferroelectric domain /? at < is replaced by the antiferro- 
electric one at > 0. These two large domains are separated by a relatively narrow 
domain within which the sublattice excitonic averages {Pfa, Pfd) change continu- 
ously (see Fig. 6) from the ferroelectric case (P^ — P^) to the antiferroelectric case 

Qualitatively the same picture we have observed also for the three-dimensional 
case. This is illustrated in Fig. 7, where the ground-state phase diagrams of the 
extended FKM are plotted for two different values of Coulomb interaction {U — 2 and 
U — A). These results indicate that ferroelectricity remains robust with respect to 
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the increasing dimension of the system, what should be important for an apphcation 
of HF solutions on a description of real three dimensional systems. 

In conclusion, we have calculated the ground-state phase diagram of the spinless 
FKM with /-/ hopping in the HF approximation with the CDW instability We have 
found that the HF solutions with the CDW instability reproduce perfectly the two- 
dimensional intermediate coupling phase diagram of the extended FKM calculated 
by CPMC method. Using this fact we have extended our HF study on cases that 
have been not described by CPMC and namely, the case of small values of /-electron 
hopping integrals, the case of weak Coulomb interactions and the three-dimensional 
case. We have found that the ferroelectric ground state with the spontaneous polar- 
ization remains stable in all examined cases. 

This work was supported by the Slovak Grant Agency for Science under grant No. 
2/7057/27 and the Slovak APVV Grant Agency under Grant LPP-0047-06. I would 
also like to acknowledge H. Cencarikova for a technical help during the preparation 
of manuscript. 

Note added. - Aher submitting this work we came to know about the work of Schneider 
and CzychoU [15] who studied the extended FKM in the limit of infinite dimensions 

and obtained results similar to ours. 
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Figure Caption 



Fig. 1. Dependence of the HF parameters nf,5f,nd,Sd, A. and Ap on the /- 
level energy Ef calculated (with step AEj = 0.005) for three different values of tf 
{tf — —0.2, —0.5, —0.8) and U — 2. Insets show the tf — —0.5 case at much higher 
resolution (the numerical data have been obtained with step AEf — 0.00005). The 
case of tf = —0.8 is analogous to tf = —0.5. 

Fig. 2. The HF (•) and CPMC [9] (□) phase diagram of the two dimensional 
FKM with /-/ hopping obtained for [/ — 2. 

Fig. 3. Dependence of the HF parameters A, Ap and 5d on the /-level energy 
Ef calculated for different values of tf (t/ = 0, -0.01, -0.02, -0.05) and [7 = 2. 

Fig. 4. The complete HF phase diagram of the two-dimensional extended FKM 
in the intermediate {U — 2) and weak coupling {U — 1) regime. 

Fig. 5. Dependence of the HF parameter A on the Coulomb interaction U calcu- 
lated for different values of tf and Uf = 1/4. 

Fig. 6. Dependence of the excitonic expectation value Pfd = {ft^i) on tf calcu- 
lated for Ef = 0.7 and U = 2. 

Fig. 7. The complete HF phase diagram of the three-dimensional extended FKM 
calculated foicU — 2 and U — A. 
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